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We analytically compute the large-deviation probability of a diagonal 

H^H 



U ' 



matrix element of two cases of random matrices, namely /3 = 

and 7 = \1n + pH^H] , where H is a M x iV complex Gaussian matrix 
with independent entries and M > N. These diagonal entries are related 
to the "signal to interference and noise ratio" (SINR) in multi-antenna 
communications. They depend not only on the eigenvalues but also on 
the corresponding eigenfunction weights, which we are able to evaluate on 
average constrained on the value of the SINR. We also show that beyond a 
lower and upper critical value of /3, 7, the maximum and minimum eigen- 
values, respectively, detach from the bulk. Responsible for this detachment 
is the fact that the corresponding eigenvalue weight becomes macroscopic 
(i.e. 0(1)), and hence exerts a strong repulsion to the eigenvalue. 

1. Introduction 

Random matrix theory has recently seen a flurry of applications in com- 
munications. Here, the random matrix under study may be the matrix 
of random channel amplitudes between transmitting and receiving multi- 
antenna arrays, [HE] or an array of pseudo-random code vectors used in 
a multi-user communications setting in order to scramble the signals from 
other users [3l[4J. 

One metric to characterize the performance of the communications is 
the mutual information, which gives the ultimate number of bits per chan- 
nel use that can be transmitted without error. The ergodic mean and the 
fluctuations of this quantity has been analyzed under a wide range of as- 
sumptions regarding the channel statistics [2l[5l[6]. For slowly time- varying 
channels a better metric for the performance is the so-called outage capac- 
ity which provides an achievable transmission rate given a probability that 
this rate will not be achieved from the underlying fading channel [7]. As a 
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result, a number of works showed that for large antenna numbers the fading 
statistics become Gaussian [8,9j. More recently, the tails of the distribution 
were also calculated using the Coulomb Gas approach |10j . 

To obtain the above full advantages from multiple antennas, it is neces- 
sary to have an optimal receiver structure, which however is quite complex 
to implement in real systems. Instead, low complexity albeit suboptimal 
linear receivers offer practical alternative. 

Such receivers include the so-called MMSE (minimum mean square er- 
ror) and the zero-forcing (ZF) receivers. The information throughput per- 
formance depends on the ability of the linear receiver structure to mitigate 
interference. One very useful method to quantify the performance is through 
the asymptotic analysis of the signal to interference and noise ratio (SINR) 
for the MMSE receiver in the limit of large antenna numbers using tools 
from random matrix theory. 

As in the case of the mutual information, when the channel is slowly 
varying, it is important to evaluate the full probability distribution of the 
SINR to obtain the probability of outage for a given target SINR value. 
This is important when the number of antennas is not too large, in which 
case the fluctuations play an important role. 

In a seminal work [11] the authors proved the asymptotic normality of 
the SINR for the MMSE and ZF receivers when all transmitters have equal 
power. More recently, [12|I13] showed the normality of the MMSE SINR. 
Unfortunately and in contrast to the total mutual information, the Gaussian 
approximation for the SINR behaves badly unless the number of channels 
is quite large. As a result, inspired by the fact that the SINR for the equal 
power MIMO ZF receiver has a Gamma distribution [111114) . several works 
were devoted in approximating the SINR statistics with other distributions, 
notably the Gamma and generalized Gamma probability densities [15H18j . 
by matching the first three moments. Nevertheless, this methodology, al- 
though perhaps providing good agreement under certain conditions, is ad- 
hoc and does not offer any intuition on the SINR statistics. Finally, it 
should be pointed out that the exact distribution of the SINR has been cal- 
culated recently in terms of ratios of determinants [19j. Nevertheless, such 
an analysis is quite tedious and does not provide any intuition about the 
result. 

In this paper, we take a different approach. Instead of trying to prove 
Gaussian behavior close to the peak of the distribution of SINR, we apply 
the Coulomb Gas methodology, which allows us to calculate the distribu- 
tion of the SINR arbitrarily far from its most probable, ergodic value. The 
Coulomb gas model was introduced originally by [20] and more recently has 
seen numerous applications [IU 1[^TH25] treats each eigenvalue of a random 
matrix as a point charge in the presence of an external potential while re- 
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pelling the others. To apply this model, we rely on the fact that the SINR 
can be written in terms of a diagonal matrix of a random matrix and hence 
as a sum over the eigenvalues of the matrix. Nevertheless, since this sum de- 
pends not only on the eigenvalues but also the weights of the corresponding 
eigenfunctions on the matrix element, we need to generalize the Coulomb 
gas approach to take into account the effects of the fluctuating weights. It 
should also be mentioned that single matrix element distributions of related 
quantities have been evaluated elsewhere but in a different context 
and without exemplifying the interaction between eigenvalues and eigen- 
functions. 

Outline: In the next section we present the channel model and introduce 
the concepts of the SINR for the ZF and MMSE receiver. In Section [3] 
we present our analytical results, while in Section H] we demonstrate their 
validity numerically and we conclude in Section O 



2. Problem Statement 

In this section we define the channel model. We consider a wireless 
communications system with an N antenna transmitter array and an M 
antenna receiver array. It is typically assumed that M > N and the ratio is 
defined as a = M/N > 1. The M-dimensional received signal vector y can 
be written as 

y = Hx + z (1) 

where the vector x represents the transmitted signal with identically dis- 
tributed elements and variance E[xx^^] = plj^. z is the noise vector, with 
independent complex Gaussian elements ~ CJ\f{0, 1). The channel matrix 
H is assumed to have independent elements ~ CA/'(0, 1/N). 

The basic communications problem at the receiver is try to deduce x 
from y, given the knowledge of the channel H. Even though the optimal re- 
ceiver structure leads to the maximum throughput per channel use there are 
several suboptimal receivers, which are popular because they are linear in 
their implementation. The most common ones are the so-called Minimum- 
Mean-Square-Error (MMSE) and the Zero-Forcing (ZF) receivers. In both 
cases the vector y is multiplied by a matrix P in an effort to mitigate the 
noise and interference. 

2.1. MMSE Receiver 
In this case the matrix Pmmse is 

Ht (2) 
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This matrix has the property of minimizing the average square error of 
the signal in the presence of the noise z. The output signal can be then 
expressed as 

X — Pmmsey (3) 

The resulting signal-to-interference-and-noise ratio (SINR) for each signal 
stream Xi for i = 1 , . . . , is given by [3] 



1 



(I^ + pHtH)-^ 



1 



(4) 



It will turn out to be convenient to parameterize this quantity by Zi = ■yk/p 
and re-write the above equation as 



1 

1 + pzi 



N 



i=i 



+ pXj 



(5) 
(6) 



In the second line we have expressed the i-th diagonal element of the matrix 
in the RHS of ^ in terms of the eigenvalues xj and the matrix elements 
of the unitary matrix U, which diagonalizes H^H. Since the elements of H 
are ~ CM{0,1/N), U is Haar-unitary matrix. As a result, the quantities 
\uji\'^ for fixed i and j = 1,. . . ,N are uniformly distributed in (0, 1) with 
the constraint 



N 

El 



(7) 



2.2. ZF Receiver 

Similarly the SINR of the zero-forcing (ZF) receiver can be obtained. In 
this case the projector matrix P^j is simply the pseudo-inverse of the matrix 
H (H"*"), which if course exists with probability one only when M > N . As 
a result, the output vector is 

X = P,/y = X + H+z (8) 

We see that multiplication with on y kills all self-interference of signals, 
since all terms involving signals Xg, with q ^ i are forced to zero (hence 
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the name "zero- forcing" ) . Of course this comes at the cost of increasing the 
noise. The corresponding SINR can be written as 



P_ 



1 

Zi 
N 

E 



HtH 



(9) 

(10) 



with the quantities uji and xj defined as above. 

It is worth pointing out that the SINR of both MMSE and ZF cases 
above may be written as a sum over a function of eigenvalues weighted by 
the corresponding eigenvector weight. 



N 

E 



u 



'j'li 



S[Xj 
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Also, it is important to mention that in the hmit of large p, Zmmse coincides 
with Zzf, i.e. z^f = liiUp^oo Zmmse and thus /3{p) = />limp/_>oo 7(p')/p'- As 
a result, we will focus on the distribution of the MMSE SINR first, from 
which we will be able to derive all results for the ZF SINR by taking the 
appropriate limit. 



3. Technical Analysis 

In this section we will go through the basic steps of the calculation of the 
probability distribution function (PDF) of the normalized SINR Zi, omitting 
the index i when necessary. Keeping in mind the statistics of Uji and their 
constraint ([7]) we write the PDF of z as 

(12) 

where the expectation is over the vector x of the H^H eigenvalues and 
the random vector t with distribution identical to the quantities A^|?Xjjp, 
for fixed i and j = 1, . . . ,N. We have used the compact notation s{x) to 
indicate both ZF and MMSE cases above. For simplicity we will omit the 
dependence of s{z) on z as well as the overall proportionality factor unless 
explicitly mentioned. We may explicitly integrate over t by first expressing 
the above (5-function as well as the constraint 

N 

^t^ = N (13) 
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x,t 
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S Ns{z) - ^s(a;j)tj 
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as Fourier integrals. As a result we obtain 



F{z) oc dk dXE- 



N „Ar 



OC 



I"! 



,N{ks+X) 



n 

N 

n 



1 _ g-Af(A:s(xj)+A) 

+ A) 



(14) 



Note that that although the integral over k and A is along the imagi- 
nary line, the saddle point will lie on the real axis and hence we omit the 
imaginary i for simplicity. The expectation over x is performed with the 
eigenvalue distribution P(x) given by 



N 



P(x) OC A(x)2 H 



^M-N^-Nxj _ g-Ar^F(x) 



(15) 



with A(x) = ni>j(^i ~ ^j) Vandermonde determinant and the second 
equation being the definition of F{x). When is large, the eigenvalues 
of H^^H form a tight density which can be represented as a density p{x) 
corresponding to the quantity 



j 

As a result we may rewrite (|14|) as 



(16) 



IZ OC 



(17) 



where J Dp represents a path integral over non- negative, normalized p{x), 
F\p\ is the energy functional associated with the probability distribution of 
eigenvalues [TU1[251[27] and Eq\p] is the functional obtained from the 
exponent in 



F[p] 
Eo\p] 



dxp{x) I X — [a — 1) ln(x) — / dx'p{x') In |x — x'\ I (18) 



-{ks + \) 



dxp{x) In 



1 _ f,-N{ks{x)+\) 

ks{x) + A 



(19) 



It is crucial to point out that the functional F[p] in is multiplied by A^^ 
while Eq\p\ is only multiplied by A^. Hence the fluctuations of F\p\ will be 
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far smaller than those of i?o[p]- As a result, to leading order in we may 
first find the optimal distribution that minimizes F\p\. This distribution is 
the celebrated Marcenko-Pastur distribution given by 



Po(x) = (20) 

where the limits of the support are 

a,h={^/^±lf (21) 

Subsequently, using this Pq{x) we may find the optimal values of k and 
A that minimize £'o[j'o]- This two-tiered approach works, as mentioned 
before, because, for large N , the eigenvalue distribution has much smaller 
fluctuations compared to the fluctuations of the unitary matrix elements 
\uji\^ and k, A. 

We next analyze the above equations in two separate regimes, depending 
on whether the quantity A + ks{x) is positive. When it is, the exponential 
factor inside the logarithm of (|19p is negligible, and we may therefore omit 
it. This corresponds to the situation when all weights of the eigenvalues 
are of similar size, i.e. \uji\^ = 0{1/N), or equivalently t{xj) = 0(1). 
The analysis of this region will be discussed next in Section 13.11 When 
A + ks{x) < 0, we need to take the exponential explicitly into account, 
which we will do in Section 13.21 In this case, as we shall see, the average 
weight of one eigenvalue becomes macroscopic, i.e. t{x) = 0{N). 

3.1. Region with A + ks{x) > 

In this case, the exponential inside the logarithm of (|19p is exponentially 
small in N and therefore may be neglected. As mentioned above, all typical 
values of tj are of order unity, i.e. tj = 0(1), with their sum fixed to 
([13]). In fact, the integral over p{x) in (HHD represents the entropy of the 
random variables t for given p{x) and the mean value of the tj = A^|njjp is 
equal to 

The saddle point equations for A and k are obtained by differentiating 
EqIpq] with respect to A and k, respectively, and setting the derivative to 
zero 

dx-P^ = 1 (23) 



a A + ks{x) 
Po{x)s{x 
A + k.s{x 
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By identifying l/{ks{xj) + A) as the average value of tj, we immediately 
see that the first equation is nothing else but the normalization condition 
(fT3|) . Similarly, the second equation simply states that — 
i.e. imposes the <5-function constraint in (jl2p . We note that combining the 
two equations we get the identity 

X + ks{z) = l (25) 

This, together with e.g. (|24p will provide us with the optimal values of k, A 
to plug into (jl9p and thus evaluate the leading term in the exponent of the 
probability distribution of s (resp. z) by evaluating it at the saddle point. 

We start by making the following convenient change of variables from k 
to c through 

which for the MMSE case becomes k = — A(l + pc). Once this variable is 
determined, the other, e.g. A can be obtained from it through (j25p 



A= (27) 

s{c) - s{z) 

Plugging ()26p into ()24p we get, after some re-arrangements, 

= -±- (28) 
X — c z — c 

where we have also used the fact that s{z) = 1/(1 + pz). It is interesting to 
point out that this equation is independent of p and holds also for the ZF 
receiver. Also it represents a balance of forces for a (yet fictitious) charge 
located at c: from one side we have the repulsion of the Coulomb sea, while 
from the other there is another (fictitious) charge located at the position 
dictated by the normalized SINR z. Before proceeding to integrate the LHS 
we note that, in order to get a convergent answer, c has to take values 
outside the support of p{x), i.e. c ^ (a, &). We then have 



_ sgn(l + a - c)^/{b - c){a - c) + c + a - 1 
^ — 2 ^ ' 

We see that the region of z for which the above equation has solutions is 
— a| < ^/a. This corresponds to values of c in the regions — oo < c < a 
(for a — -y/a < z < a) and b < c < +00 (for a < z < a + y/a). Solving the 
above equation for c gives us 

c(z) = z(l + ^— ) (30) 
\ z — a J 
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The values of c{z) for which the solutions above break down are c = a 
{z(a) = a — ^/a) and c = b {z{b) = a + ^/a). 

We may now calculate the exponent of the PDF for Iz — a| < ^/a. To do 
so, we simply need to plug in the above values of k and A (obtained directly 
from c) into (jl9p and calculate the corresponding integrals. As discussed 
before p^[22| [27] the value of F\po] does not depend on z and is therefore a 
constant (we have also omitted the dependence of Eq on po). 

-1 



^o(^) 



In 



p-^+z 



+ 



c + p 



(31) 



+ 



{sgn{l + a- c)^/{b-c){a-c) - ^/ {p-^ + a){p-^ + b)'^ 



+ (a + l)ln 



V^|6-c| + ^ 



Vp~i + b + Vp-i + 



[a 



Din 



\/a\b — c\ + \/b\a 



^a{p-^ +b) + ^b{p-^ + 



Plugging in the dependence of c{z) we obtain the following simplified for- 
mula. 

^-1 



Eo{z) 



P 



1 



a 



z — alnz + ln{p + z) + 
y/{p-^+a){p-^+b) ^^^^ ^ (a + l)lna 



(32) 



(a + 1) In 



+ a + + b 



+ (a-l)ln y/b{p-^ + a) + y^a{p-^ + b) 

This exponent has two interesting properties. First, it is a maximum at the 
ergodic value of z, which corresponds to the ergodic average of the SINR. 
This corresponds to k = 0, and, following (|25p also A = 1. As a result, (|24p 
equates s{z) to its ergodic average over the Marcenko-Pastur distribution, 
Epp[s(x)]. As we shall see, this corresponds to the peak of the Gaussian 
distribution of the SINR distribution. The MMSE SINR is the given by 



Terg — P^erg^mmse 



V^(l - (a - l)p)2 + 4ap + p(q - 1) - 1 



(33) 



This can be seen by directly maximizing ()32p over z. Indeed, expanding 



(|32]) close to Zerg wc find 



i^o(^) 



■'erg j 



2v. 



(34) 
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where Verg is the variance of the MMSE SINR given by [28j 

_ (a - 1)a/(1 - (a - + 4ap + p{a - 1) + a + 1 



2^(l_(a-l)p)2 + 4a/> 



(35) 



The corresponding values of /3erg and v^rg for the case of the ZF receiver 
can be obtained by taking the hmit p ^ oo but keeping z fixed, and then 
multiplying with p to get the SINR fS^rg 

/3erg = pz^rg^zf = P{a - I) (36) 
Verg^zf = a - 1 (37) 



We see that when a = 1 the Gaussian approximation breaks down |1HI28]. 

Keeping only the dependence on z in ([32]) . we see that, to leading expo- 
nential order. 



z 



M 



(-) « i^p~i + ,)N ^~''' (38) 

This formula is remarkable for two reasons. First, it is surprising that 
this simple formula peaks at the ergodic value of Zerg in (f33|) . Second, it 
settles a year-old conjecture, that the distribution of MMSE SINR should 
be (approximately) a Gamma distribution. Several papers in the literature 
|15H17] tried to fit the distribution to the Gamma distribution by fitting 
their moments. We see that this asymptotic form illustrates that although 
simple in form, it is not a Gamma distribution. 

By letting p — )• oo we can recover the distribution of z^f 

P^/(z) oc z^^'-"e-^' (39) 
P,;(/3) oc pM-N^-Np/p 

It turns out that this result is in fact exact |1HI14). 

Finally, we can also evaluate the average weight of each eigenvalue x G 
[a, 6] constrained on the value of z. Using ([2^ . ([77]) and ([50]) we obtain 

/ , M 1 S{z) — sic) , , , 

^[t{x\z)] = — - = -^^ ^ 40 

^ ^ ' A + ks{x) s{x) - s{c) ^ ' 

for |z — a| < -y/a, which is valid for both ZF and MMSE. This result, can 
also be obtained by noting that the distribution of t{x\z) is exponential 
~ exp[— (A + ks{x))t\^ as seen in (HH). 
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3.2. Region with X + ks{x) < 

Before moving on, it is worth pointing out that the above behavior is 
bound to break down at some value of z. Indeed, assuming A + ks{x) > 
and using ([23]) . (pl|) and the fact that s{x) is a decreasing function of x, we 
get the fohowing inequahty 

s{z) 

s{z) 

Therefore, for z < a and z > b, the assumption A + ks{x) > and the 
resulting equations ([23|) and (pl]l have to break down. 

To see how, we need to analyze the situation outside the region — a| < 
y/a. We thus need to consider the situation when for some eigenvalue(s) 
the exponent in ()19p becomes positive, i.e. when A + ks{x) < 0. For the 
section of the support of p{x) where this occurs, the exponent in (jl9p will 
be positive, so we will need to include an additional term in -Eobl) namely 

NEolp] ^ -{ks{z) + X)N + N j dxp{x)\n\ks{x) + \\ (42) 

J a 

+ dxp{x){ks{x) + X) 

Jn 

where TZ is the region of the support of p{x) (possibly including only a finite 
number of eigenvalues) with X + ks{x) < 0. This extra term is an additional 
potential of strength N'^ks{x) exerting a force on the charge density p{x). 
Since it is 0{N'^) we can no longer assume it is small and we have to 
take it into account explicitly together with N'^F[p\ in the determination 
of the optimal p{x). First, we need to estimate whether the number of 
eigenvalues affected is finite or scales with N . To answer this we start by 
observing that for these eigenvalues the corresponding typical value of tj 
becomes of 0{N). Due tot the constraint Ylij'^j — ^ (fT3|) . there can be 
at most a finite number of such eigenvalues with corresponding tj = 0{N). 
We will initially assume that it is only one such eigenvalue and later on 
show that this is consistent. As a result of these considerations we need to 
treat this eigenvalue separately from the others. Therefore, we separate the 
continuous part of the eigenvalue density and express is as 



^ p{x)s{x) 



X + ks{x] 
^ ^ p{x)s{x) 



X + ks(x] 



< s{a) 
>s{b) 



(41) 
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and denote the position of the eigenvalue (which may be the largest or 
smallest depending on whether we are analyzing the case z > a + \foL or 
z < a — i/a, respectively) by y. The exponent in p!7|) can be expressed as 

-{N -\fF\q\-NEM (44) 

where F\q\ is the same energy functional as in (jl7p . As a result, the optimal 
q{x) is still the Marcenko-Pastur distribution (|2U|) . Thus £'+[po] is given by 

y] = ?/ — (a — 1) Iny — 2 / dxpQ{x)\a.\x — y\ — (ksi^z^ ^ (45) 

J a 

b fb 

dxpoix) ln{ksix) + A) - / dxpoix) In 1 - e-^('=^(^)+^) 

J a 

+ ^ (in |A + ks{y)\ - In ^e-^i''^iy)+^) _ i 

This need of explicitly splitting one eigenvalue from the bulk and treating it 
in a special way has appeared also in the context of bipartite entanglement 
|25y29j. Since only ks{x) + A > for x ^ y the last term in the second line 
above will only contribute subleading terms and therefore may be neglected. 
We now need to find the saddle point jointly for y, A and k. 

^ ~ \ + ks{x) ^ N{\ + ks{y)) ^ e-^^^+^<y)) - 1 ^ ' 

J a A + ks{x) ^ N{\ + ks{y)) ^ e-m^+f'sm _ i > 

'dxP^ + ^ (48) 
y -X y 

kps{yf kps{yfe-^^^+^<y^^ 
^ N{\ + ks{y)) ^ e-^i^+f^'iy)) - 1 

From the first equation we conclude that the values of A + ks{y) = 0(1/N). 
Otherwise, if A + ks{y) = 0(1), the integral in the right-hand-side of (|46p 
would have to vanish, which is inconsistent with the fact that \ + ks{x) > 0. 
Thus, setting X + ks{y) = —w/N for w still unknown, we find that to leading 
order, 

k = . . \ . (49) 



and tt; is a solution of the equation 

.^)is{z) 

s{x) — s{y) w 



^,«(£M£>^.i l_ (50) 



13 



Putting this together, we finally get the equation for y: 

p{x)dx _ ^ g - 1 l_ ^^^^ 

Ja y-x y y-z 

This is last equation is both surprising and intuitive. It is firstly surprising 
that the value of y does not depend on the specific form of the SINR function 
s{z) but only on z itself, the normalized SINR. Second, it tells us that the 
position of this eigenvalue is determined by external forces exerted on the 
other eigenvalues and, in addition, it feels the repulsion of a unit charge 
located in the position z. Also it is interesting to note that it feels only half 
the repulsion from the Marcenko-Pastur continuous eigenvalue density. The 
other half has been screened away due to the interaction of this eigenvalue 
with the matrix elements of the diagonalizing unitary matrix (through A and 
k). It should be pointed out that ([5T]) is valid for both cases z < a — ^/a 
and z > a + -y/a, with y < a and y > b, respectively. 

In the above analysis we have assumed there is only one eigenvalue that 
detaches from the bulk. Let us assume there were r > 1 such eigenvalues. 
In that case, they would all have to stick together satisfying A + ks{yj) = 
—Wj/N. Otherwise, if say only one yi satisfied this relation, all others with 
yj — yi = 0(1), for j = 2, . . . ,r, would necessarily have A + ks{yj) > 0, 
which would not be sufficient to provide the "kick" to get out of the bulk. 
However, on the other hand, if these r eigenvalues are within 0{1/N) from 
each other, their repulsion l/{yi — yj) will be large {0{N)) and hence would 
dominate (05]). As a result, only one eigenvalue can be detached from the 
bulk. 

We may now integrate the LHS of (|5ip to get 



^/{a-y){b-y) + y-a + l ^ ^ _ a-l l_ 

2y y y- z 

Solving for y gives 

y = z(l + ^— ) (53) 
V z — a I 



\z — a\ > va 



which is identical with (j30p . although obtained through a completely dif- 
ferent method, and with z here taking different values. One way to jointly 
interpret c and y is that they correspond to the location of a "state" , which 
when located outside the continuum of nearby states forms a bound state 
(y), while when it enters the continuum, it becomes a "resonance". 



14 



We may now plug in the above results into ()45p to obtain the exponent 
of the PDF in this region of z. The final result is 



E^{z) 



In 



V\ 



+ 



-1 



(a - l)lny 



+ 



p-^ + z_ 

^ {sgn{y - a)^/{y-h){y-a) - ^/ {p-^ + a){p-^ + h)) 

+ In 



(54) 



- (a + 1) (\u\^/\y-h\ + ^\y-a 
+ (a — 1) ^In \/ a\h — c| + \/h\a — 
+ In v^a(p-i + 6) + v^6(p"i + a) ) 

Inserting ()53p into ()54p we find that, up to a constant, Ej^{z) is identical 
to Eq{z) in ()32p . thus extending the validity of the latter for all values of 
^ > 0. 

Finally, we evaluate the average eigenvalue weights constrained on the 
value of SINR (or equivalently to z), which are the analogues of (00]). Using 
19]) we find the identical expression as (|40p 



¥.[t{x\z)] 



1 
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In addition, we may calculate the mean weight of the detached eigenvalue 
y. As with the other weights, its distribution is exponential ~ exp[— i(A + 
^^(y))]) (mi)- We thus find that in this case 
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z{{l + pa){z — a) + pa) 
(1 + pz){z - afiz -a + 1) 



for |z — a| > -^/a, where w = —N{\ + ks{y)) appears in (fSUp . We see that 
in the limit \z - a\ = ^a, Wyi^)] = (i-e- = C>(l/iV)). We plot 



4. Numerical Simulations 

To test the applicability of this approach, we have performed Monte- 
Carlo simulations and have compared our large deviations Coulomb Gas 
(CG) approach with Monte Carlo (MC) simulations and the Gaussian ap- 
proximation. In Fig. 2 we plot the normalized probability density of the 
SINR for small p = 1, while in Fig. 3we plot it for larger p = 10. In both 
cases we see good agreement. 
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(a) Weight for ZF (b) Weight of x^^„ 

Fig. 1. (a) Average weight of the eigenfunction j on the i element, E[iV|uji p] for the 
ZF function of the corresponding eigenvalue Xj , constrained on different 

values of z. The five curves include to the lower critical {z = a — \/a), upper 
critical (z = a -\- ^/a) and an intermediate (z ^ a) value of z. The remaining two 
curves have z below the lower critical value z — a — ^/a — 0.5 and above the upper 
critical value z — 2a. In the first two we clearly see the divergence at the lower and 
higher edges of the spectrum, corresponding to the fact that beyond these values 
the weight of the edges becomes macroscopic 0{N). (b) Average weight of the 
minimum eigenvalue for < z < a — \fa. Note that now we plot the macroscopic 
occupation of the eigenvalue, i.e. E[|umi„,fep]. We plot the case of MMSE for 
various p as well as the case of ZF. 



5. Conclusion 

In this paper we have used a large deviation approach to calculate the 
probability density of the "signal to interference and noise ratio" (SINR) 
for multi-antenna arrays for two popular receiving algorithms, namely the 
MMSE and the ZF algorithms. The approach is formally valid for large N 
antenna numbers, but is not restricted to the behavior close to the peak of 
the distribution, which has been shown to be asymptotically Gaussian when 
the number of antennas is very large. Instead we calculate the probability 
of the SINR being arbitrarily away from its ergodic peak. Surprisingly, the 
leading term of the exponent of the distribution is very simple, and the 
distribution is neither Gamma, nor Beta and certainly not Gaussian. In the 
ZF case, we recover the known chi-square result. We also test the MMSE re- 
sults numerically and find good agreement even for relatively small antenna 
arrays. From a technical point of view, since the SINR of the two algorithms 
are related to the diagonal matrix elements of the matrices 
and (H^^H) ^, the task is to find the distribution of a single diagonal matrix 
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(a) a = 1 (b) a = 2 

Fig. 2. PDF of MMSE SINK for M = 6, p = 1. The agreement of the Coulomb 
Gas (CG) curve with Monte-Carlo (MC) simulations is good, even for such small 
matrices, especially compared to the Gaussian approximation. Denoted with circles 
are the values of z at which the "inner" and "outer" solutions match and we see 
no discontinuity in the numerics. 




element. The methodology we applied is based on the so-called Coulomb 
Gas model, in which each eigenvalue can be seen as a point charge interact- 
ing with an external potential and repelling each other. In this particular 
case however, the eigenvalues interact not only with each other but also with 
the weights of their corresponding eigenfunctions in the particular matrix 
element. As a byproduct of our analysis we are able to calculate the aver- 
age weight of each eigenvalue in the particular matrix element, constrained 
on the value of the SINR or the matrix element. The interaction between 
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the eigenfunction weights and the corresponding eigenvalues can be quite 
strong and as a result, below and above critical values of z the lowest and 
largest eigenvalues detach from the bulk. Nevertheless, it seems that there 
is no discontinuity involved in this detachment, at least to leading order. In 
hindsight, this is not surprising. A given diagonal matrix element depends 
on a number of 0{N) random variables of the matrix, which has 0{N'^) 
random variables. In our approach we have expressed this diagonal matrix 
in terms of the eigenvalues, which depend on the whole matrix. Somehow, 
we expect that the interaction with the eigenvalue weights will "wash" out 
this dependence from the full matrix. 
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